Nonadditive entropy reconciles the area law in quantum systems 
with classical thermodynamics 
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The Boltzmann-Gibbs-von Neumann entropy of a large part (of linear size L) of some (much 
larger) d-dimensional quantum systems follows the so-called area law (as for black holes), i.e., it is 
proportional to L d_1 . Here we show, for d = 1,2, that the (nonadditive) entropy S q satisfies, for a 
special value of q ^ 1, the classical thermodynamical prescription for the entropy to be extensive, 
i.e., S q oc L d . Therefore, we reconcile with classical thermodynamics the area law widespread in 
quantum systems. Recently, a similar behavior was exhibited, by M. Gell-Mann, Y. Sato and one 
of us (C.T.), in mathematical models with scale-invariant correlations. Finally, we find that the 
system critical features are marked by a maximum of the special entropic index q. 
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The aim of statistical mechanics is to establish a direct 
link between mechanical microscopic laws and classical 
thermodynamics. The most famous classical theory in 
this field was developed by Boltzmann and Gibbs (BG) 
and it is considered as one of the cornerstones of con- 
temporary physics. The connection between the micro- 
and macro- world is usually described by the so-called BG 
entropy, 
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where k is a positive constant, W is the number of micro- 
scopic states, and {pi\i=x,...,W is a normalized probability 
distribution. The subtle concept of entropy lays the foun- 
dation of classical thermodynamics. The BG entropy is 
additive, i.e. Sbg{A b ) = Sbg(A) + S B g{B), where A 
and B are two probabilistically independent subsystems. 
One of the crucial properties of the entropy in the con- 
text of classical thermodynamics is extensivity, namely 
proportionality with the number of elements of the sys- 
tem, when this number is large. The BG entropy satisfies 
this prescription if the subsystems are statistically (quasi- 
independent, or typically if the correlations within the 
system are generically local. In such cases the system is 
called extensive. 

In general, however, the situation is not always of this 
type and correlations may be far from negligible at all 
scales. In such cases the BG entropy (of the entire sys- 
tem or large part of it) may be nonextensive. Nonethe- 
less, for an important class of such systems, an entropy 
exists that is extensive in terms of the microscopic prob- 
abilities [U . The additive BG entropy can be generalized 
into the nonadditive g-entropy 0, [ 
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This is the basis of the so-called nonextensive statistical 
mechanics 0,0] (see H| for a regularly updated bibliog- 
raphy), which generalizes the BG theory. In the context 



of cybernetics and information theory, the same type of 
entropic form has been advanced in [6|. 

Additivity (for two probabilistically independent sub- 
systems A and B) is generalized by the following pseudo- 
additivity: S q (A,B)/k = S q {A)/k + S q {B)/k + (1 - 
q)S q (A)S q (B)/k 2 . For subsystems that have special 
probability correlations, extensivity is not valid for Sbg, 
but may occur for S q with a particular value of the index 
q =/= 1, called the g-entropic index q en t- Such systems are 
sometimes referred to as nonextensive The phys- 

ical relevance of the g-entropy is related also to the re- 
cent generalization of the central limit theorem (q-CLT) , 
where g-Gaussian distributions take the place of the usual 
Gaussians when considering systems with strong correla- 
tions Q ■ Much experimental evidence of predictions of 
nonextensive statistical mechanics is available: see, for 
instance, [1, Q . Let us emphasize the difference between 
additivity and extensivity for the entropy. Additivity de- 
pends only on the mathematical definition of the entropy; 
therefore, Si is additive, whereas 5^ (q ^ 1) is nonaddi- 
tive. Extensivity is more subtle, since it also depends on 
the specific system, as we will show in this paper. 

A physical system may exhibit genuine quantum as- 
pects. In particular, quantum correlations, quantified 
by entanglement, can be present. The classical proba- 
bility concepts are replaced by the density matrix op- 
erator p, in a more general probability amplitude con- 
text. The quantum counterpart of the BG entropy, which 
is called the von Neumann entropy, is thus given by 
Si(p) = — fcTr/Sln/5, while the classical g-entropy is re- 
placed by: 



1 - Tr p q 
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The pseudo-additivity is now given by 



(3) 



S q {pi®h) = Sg(pi) S q (p 2 ) , 5 g (pi) S q (p 2 ) 

h h h ^ ' h h 

from now on k=l. 
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In the following we will analyze fermionic and bosonic 
quantum systems in which strong non-classical correla- 
tions exist between their components. The appearance of 
long-range correlations in the ground state of a quantum 
many-body system, undergoing, for instance, a quantum 
phase transition at zero temperature, is deeply related to 
entanglement [To| . Quantum spin chains, composed of 
a set of localized spins coupled through exchange inter- 
action in an external transverse magnetic field, capture 
the essence of these intriguing phenomena and have been 
extensively studied p], El El El EE El ■ Th e degree 
of entanglement between a block of L contiguous spins 
and the rest of the chain in its ground state is measured 
by the von Neumann entropy of the block. For large 
block size, it typically saturates off criticality, whereas it 
is logarithmically unbounded at the critical point. 

Here we show that the nonadditive entropy d HJ 
S q (pL) = k 1 g T 'f ^ , Pl being the density matrix of a 
block of L spins of the ground state of quantum spin 
chains in the neighborhood of a quantum phase transi- 
tion, is extensive [i.e., S q (pl) oc L for I » 1] for spe- 
cial values of q < 1. The additive von Neumann en- 
tropy S\(pl) = — fcTrpilnpL is nonextensive [indeed, 
linii^oo Si(pl)/L = 0] in all considered cases; the same 
happens with the additive Renyi entropy. 

A similar behavior can be observed for another impor- 
tant class of quantum systems, harmonic lattice Hamilto- 
nians, i.e., bosons interacting through a quadratic Hamil- 
tonian in dimension d = 2. These systems are discrete 
versions of the free scalar Klein-Gordon field. The degree 
of entanglement between a square block of L 2 oscillators 
and its exterior, as measured by the additive von Neu- 
mann entropy of the square block, Si (pj,2 ) is linear in 
L (area law), i.e. it is nonextensive |l7l. Ila. Il9l. I20L l2l| . 
Here we show that the nonadditive entropy S^p^) is ex- 
tensive [i.e., S q (pL2) oc L 2 for L> 1] for special values 
of q < 1. 

Therefore, we present here two physical realizations, in 
many-body Hamiltonian systems, of the abstract mathe- 
matical examples recently exhibited by M. Gell-Mann, Y. 
Sato and one of us (C.T.) in Ref. [l|, that, for anomalous 
values of q, the nonadditive entropy S q , can be extensive, 
as expected from the Clausius thermodynamical require- 
ment for the entropy. In other terms, we show explicitly 
the first two physical examples of the fact that the non- 
additive entropy reconciles the area law (e.g., typical of 
black holes) with classical thermodynamics. In addition 
to the considerable advantage of S q enabling thus the 
use of all standard thermodynamical relations, this con- 
stitutes a powerful tool to detect strong nonlocal correla- 
tions in quan tum many-body systems, by a nonadditive 
measure [H, [H U ||E El, Hi H|. Indeed, this en- 
tropic index q presents a maximum when the correlation 
length is divergent in the ground state of these quantum 
systems. 



I. QUANTUM SPIN CHAINS 

First of all, we focus our investigations on a one- 
dimensional spin-1/2 ferromagnetic chain with an ex- 
change (local) coupling and subjected to an external 
transverse magnetic field, i.e., the quantum XY model. 
The Hamiltonian is given by: 

N-i 

« = - £ [a + iPpj+i + (i - 7)*M+i + 2 M] ( 4 ) 

where &"(a = x, y, z) are the Pauli matrices of the j'-th 
spin, N is the number of spins of the chain, and 7 and A 
characterize, respectively, the strength of the anisotropy 
parameter and of a trasverse magnetic field along the z 
direction. This model for < I7I < 1 belongs to the Ising 
universality class and it actually reduces to the quantum 
Ising chain for | -y | = 1 . At T = 0, this system undergoes 
a quantum phase transition at the critical point |A C | = 1 
in the thermodynamic limit N — > 00. For 7 = it is the 
isotropic XX model, which is critical for |A| < 1 [Toj . 

Entanglement in the neighborhood of the quantum 
phase transition has been recently widely investigated 
[H El El El, EE El- In particular it has been 
shown that one-site and two-site entanglement between 
nearest or next-to-nearest spins display a peak near or 
at the critical point [ll|, E3|- Moreover, the entangle- 
ment between a block of L contiguous spins and the 
rest of the chain in the ground state, quantified by the 
von Neumann entropy, presents a logarithmic divergence 
with L at criticality, while it saturates in a non-critical 
regime [HEIEH- 

The inadequacy of the additive von Neumann entropy 
as a measure of the information content in a quan- 
tum state was pointed out in Ref. [23| , A theoreti- 
cal observation that the measure of quantum entangle- 
ment may not be additive has been discussed in Refs. 
[H HI 11 HE HI, [S3, HI . Recently, Ref. 29] suggested 
to abandon the a priori probability postulate going be- 
yond the usual BG situation. 

Here we propose to extend the definition of the von 
Neumann entropy to a wider class of entropy measures 
(see also Ref. [30|) which naturally include it, thus gen- 
eralizing the notion of the block entanglement entropy. 
The block g-entropy of a block of size L is simply defined 
as the g-entropy, Eq. ©, of the reduced density matrix 
/3l of the block, S q (pi,), when the total chain is in the 
ground state. In the following we show that, in contrast 
to the von Neumann entropy, there exists a q value for 
which S q (pL) is extensive. This value does depend on the 
critical properties of the chain and it is consistent with 
the universality hypothesis. 

The XY model in Eq. @ can be diagonalized exactly 
with a Jordan- Wigner transformation, followed by a Bo- 
goliubov rotation [l5l l3ll . [H [33| ; this allows one to an- 
alytically evaluate the spectrum of pl in the thermody- 
namic limit N — > 00. More details are shown also in Ref. 
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FIG. 1: Block g-entropy S q (ph) as a function of the block size 
L in a critical Ising chain (7 = 1, A = 1), for typical values of 
q. Only for q — q ent — 0.0828 is s q finite (i.e., S q is extensive); 
for q < q ent (q > q en t) it diverges (vanishes). 



We first analyze the anisotropic quantum XY model, 
Eq. (J4j) with 7 7^ 0, which has a critical point at A c = 1. 
The block g-entropy as a function of the block size can 
show completely different asymptotic behaviors, when 
the entropic index q is varied. In particular, here we 
are interested in a thermodynamically relevant quantity, 
namely the slope, denoted s ? , of S q versus L. It is gener- 
ically not possible to have a finite value of s\\ the entan- 
glement entropy, evaluated by the von Neumann entropy, 
either saturates or diverges logarithmically in the ther- 
modynamic limit, for respectively non-critical or criti- 
cal spin chains [lj, [13, [l~5| . The situation dramatically 
changes when the q-entropy formula in Eq. ([3]) is used: 
qualitatively it happens that, regardless the presence or 
absence of criticality, a A-dependent value of q, denoted 
q e nt, exists such that, in the range 1 <C L <C £ (£ being 
the correlation length) , s qent is finite, whereas it vanishes 
(diverges) for q > q ent (q < q en t)- We note that here the 
nonextensivity (i.e., q ^ 1) features are not due to the 
presence of say long-range interactions [H, [36| but they 
are triggered only by the fully quantum nonlocal corre- 
lations. In Fig. [1] we show, for the critical Ising model 
(7 = 1, A = 1), the behavior of the block q-entropy 
with respect to the block size: S q (pL) becomes extensive 
[i.e., < IvaiL-too S q (pL) / L < 00] at a q-entropic index 
q ent ~ 0.0828 ± 10~ 4 (with a corresponding entropic den- 
sity s gellt « 3.56 ±0.03), thus satisfying the prescriptions 
of classical thermodynamics. 

A very similar behavior is shown for the non-critical 
Ising model, as well as for critical and non-critical XY 
models with < 7 < 1. The value of q en t, for which 
Sq(pL) is asymptotically extensive, is obtained by maxi- 
mizing numerically the linear correlation coefficient r of 
SqipL), in the range 1 < £ < ^, with respect to q, 
as shown in the inset in Fig. O Let us stress that, at 



FIG. 2: The A-dependence of the index q en t in the Ising (7 = 
1, solid circle) and XY (7 = 0.75, empty circle) chains. Inset: 
determination of q ent through numerical maximization of the 
linear correlation coefficient r of S q (pL). The error bars for 
the Ising chain are obtained considering the variation of q en t 
when using the range 100 < L < 400 in the search for S q (pl) 
linear behavior. At the present numerical level, a vanishing q 
cannot be excluded off-criticality because of finite-size effects. 



precisely the critical point, £ diverges, hence L is unre- 
stricted and can run up to infinity. The index q ent de- 
pends on the distance from criticality and it increases as A 
approaches A c (Fig. \2§ . It is worth stressing that our nu- 
merical results satisfy the duality symmetry A — ► 1/A, 
investigated in Ref. [37j . For the sake of clarity, notice 
that the value of q en t may be vanishing off criticality and 
this could not be found numerically because of the pres- 
ence of finite-size effects. We have also checked other 
values of 7 for the XY model and the results are very 
similar to those presented here. This fact is consistent 
with the universality hypothesis. On one hand, the XY 
and Ising models (Ising universality class) have the same 
behavior as regards the extensivity of S q (pL). In Fig. [3] 
we report the variation of s qant with respect to A. On 
the other hand, for the isotropic XX model (7 = 0) in 
the critical region |A| < 1 we find q ent — 0.15±0.01 for 
which S q (pL) becomes extensive. 

Reference [38] enables us to analytically confirm (only 
at the critical point) our numerical results. The con- 
tinuum limit of a (l+l)-dimensional critical system is 
a conformal field theory with central charge c. In this 
quite different context, the authors re-derive the result 
Si(pl) ~ (c/3)lnL for a finite block of length L in an 
infinite critical system. To obtain the von Neumann en- 
tropy, they find an analytical expression for Trp^, namely 
Trp q L ~ L- c / & {i- 1 /i) . Here, we use this expression quite 
differently. We impose the extensivity of S q (pL) and we 
find the value of q for which —c/Q(q ent — l/q en t) = 1 , *-e., 
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FIG. 3: A-dependence of the g-entropic density s 9c „ t in the 
Ising (7 = 1, solid circle) and XY (7 = 0.75, empty circle) 
models. For A = 1, the slopes are 3.56 and 2.63, for 7 = 1 
and 7 = 0.75, respectively. Notice the considerable variation 
in the values of the slope for even slightly different values of 
A out of criticality. 



FIG. 4: q en t versus c with the g-entropy, S q (pL), being exten- 
sive, i.e., limL^oo S 2 g (pz,)/Z < 00. When c increases 

from to infinity, q en t increases from to unity (von Neu- 
mann entropy) . For the critical quantum Ising and XY models 
c — 1/2 and q en t = \/37 — 6 ~ 0.0828, while for the critical 
isotropic XX model c — 1 and q en t = VTO — 3 ~ 0.16. 



consequently, linii^oo S v /^2_ 3 (pl)/L < 00. When c in- 
creases from to infinity (see Fig. [4j , q ent increases from 
to unity (von Neumann entropy). It is well known 
that for the critical quantum Ising and XY models the 
central charge is equal to c = 1/2 (indeed they are in 
the same universality class and can be mapped to a free 
fermionic field theory). For these models, at A = 1, 
the value of q for which S q (pL) is extensive is given by 
Qent — \/37 — 6 ~ 0.0828, in perfect agreement with our 
numerical results in Fig. [2] The critical isotropic XX 
model (7 = and |A| < 1) is, instead, in another univer- 
sality class, the central charge is c = 1 (free bosonic field 
theory) and S q (pL) is extensive for q ent — \/T0— 3 ~ 0.16, 
as found also numerically. Therefore, the universal be- 
havior of the g-entropic index q ent is strictly related to 
the universal role played by the central charge in confor- 
mal field theory. Eq. ([5]) represents an additional con- 
nection between nonextensive statistical mechanical con- 
cepts and BG statistical mechanics at criticality. Sec 
Ref. f39j for another connection, where once again we 
verify that the g-entropic index typically characterizes 
universality classes. Let us note that, when the critical 
one-dimensional (ID) system is a semi-infinite chain, one 
has to replace c with c/2 in Eq. © [381 ] . 

It is worth mentioning that the Renyi entropy of a 
block of critical XX spin chains has been derived ana- 
lytically in Refs. [i^, EH]. Since the Renyi entropy is 
simply connected to the g-entropy, it is possible to re- 
derive q en t for the critical XX model also from that an- 
alytical expression. Finally, let us point out that the 
reduction of the pure ground state of the full chain (at 
T = 0) to a finite block of L spins results in a mixed 
state with quantum fluctuations. A mapping of this sub- 
system within a zero-temperature XX infinite chain to a 



finite system which is thermalized at some finite temper- 
ature has been recently exhibited [HJ , thus defining an in- 
dependent effective temperature of the block. The use of 
a non-Boltzmannian distribution (e.g., the one emerging 
within nonextensive statistical mechanics) might enable 
the definition of an effective temperature which would 
not depend on L, as physically desirable. Indeed, this 
approach has been successfully implemented for e — e + 
collision experiments (43J. 



II. 2D BOSONIC SYSTEMS 

Now we present a second physical realization of the 
extensivity of S q in a bosonic 2D system at T = 0. We 
start from a bidimensional (square lattice) system of in- 
finite coupled harmonic oscillators studied in Ref. [13], 
with Hamiltonian 

where <& x ,y, ^x.y and u>o are the coordinate, momen- 
tum and self-frequency of the oscillator at site r = 
(x,y). The system has the dispersion relation E(k) = 

) + 4sin 2 k x /2 + 4sin 2 k y /2, i.e., a gap ojq at k — 

0. Applying the canonical transformation bi — y/~^ x 

($ ; + illi) with u = \/uJo + 4 and i = l,...,L 2 , the 
Hamiltonian in Eq. ^ is mapped to the quadratic 
canonical form 



llAijOj + i (a] By a] + h.c.) 



(6) 



where are bosonic operators and i, j range from 1 to 
L 2 . In Ref. [17f the authors find an asymptotic linear be- 
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FIG. 5: Block (/-entropy S q (pL) as a function of the square 
block area L 2 in a bosonic 2D array of infinite coupled har- 
monic oscillators at T = 0, for typical values of q. Only for 
q = qent ^ 0.87, s q is finite (i.e., S q is extensive); for q < q en t 
(q > <?ent) it diverges (vanishes). Inset: determination of q e „t 
through numerical maximization of the linear correlation co- 
efficient r of S q (p L ) in the range 400 < L 2 < 1600. 



havior (i.e., area law) of the block entanglement entropy 
as a function of the linear size L (they consider square 
blocks of area L 2 ), for several u>o (limi^oo Si(pl)/L 2 = 
0), no matter how close the gap energy is to zero. 

Here we study, instead, the behavior of the block q- 
entropy of the reduced density operator of a square block 
as a function of its area L 2 , when the bosonic infinite 
two-dimensional system is in its ground state. We follow 
a similar procedure to the one used above for quantum 
spin chains (see Ref. 17[ for more details). In Fig. [5] 
we show, in the case of luq — 0.01, that S q (pL) becomes 
extensive (i.e., < hin^^oo S q (pi,)/L 2 < oo) at an in- 
dex q ent ~ 0.87 (with a corresponding entropic density 
s qent 0.011; see inset in Fig. EJ)- A very similar behav- 



ior is shown for other values of the gap. Let us point out 
that, unlike the linear behavior (for any gap energy) of 
the von Neumann entropy, now the index q ent depends on 
the gap and therefore measures the presence of a progres- 
sively divergent correlation length, as shown in Fig. [6] 



III. FINAL REMARKS 

We present two quantum many-body Hamiltonian 
physical realizations of the mathematical probabilistic 
models with scale-invariant correlations recently shown 
by M. Gell-Mann, Y. Sato and one of us (C.T.) in Ref. 
[l|, in which the nonadditive entropy S q can be applied 
successfully (i.e., satisfying the classical thermodynamic 
requirement of extensivity). In this basic manner, we 
reconcile the entropy area law characterizing many quan- 
tum systems with classical thermodynamics. In addition 
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FIG. 6: ojo-dependence of the index q ent in a bosonic 2D array 
of infinite coupled harmonic oscillators at T = 0. Inset: the 
wo-dependence of the g-entropic density s 9ent . 



to that, the present results show clearly the difference 
between additivity and extensivity for the entropy. Addi- 
tivity depends only on the mathematical features of the 
entropy, e.g, Si is additive while S q (q ^ 1) is nonad- 
ditive. Extensivity is a more subtle concept and relies 
on both the mathematical features of the entropy and 
the specific physical system. Indeed, the T = block 
entropies of the 1/2-spin d = 1 quantum system at crit- 
icality are given by S\(L) cx InL (i.e., nonextensive) , 
and S,,/ n , „a _ 3 y c (L) cx L (i.e., extensive). Moreover, the 



'[v/9H 



T = block entropies of the d = 2 bosonic system are 
given by Si(L) oc L (i.e., nonextensive), and S q (L) cx L 2 
for special values of q < 1 (i.e., extensive); for instance, 
q ent ~ 0.87 for vanishing gap energy. 

More generally, it is known (see Refs. [13, [3 and 
references therein) that, for d-dimensional bosonic sys- 
tems (e.g., a black hole (T^. I20L [2l| ) . Si follows the area 
law, i.e., Si(L) cx L d ^ 1 (nonextensive). Let us point out 
that the behavior of the block entropy Si(L) for these 
quantum systems matches well-known results in confor- 
mal field theory (as noted above), where the analogous of 
the block entropy is the so-called geometric entropy, de- 
fined in the continuum [l!| [HI, 38] . As first suggested 
by 't Hooft (1985) and later shown by Callan and Wilczek 
(1994), the geometric entropy is the first quantum correc- 
tion to a thermodynamical entropy, which reduces to the 
Bekenstein-Hawking entropy for black holes [Til, [2(| [H| • 
A direct connection between entropy and boundary area 
has been clearly suggested and numerically implemented 
in Refs. [iJJ, Ho| (compare with Refs. [13, EH and refer- 
ences therein) . Finally, the relation among entanglement 
entropy, the black hole area law and other concepts such 
as the holographic bound is still an open problem 44j . In 
this context, it is interesting to note that a logarithmic 
behavior for d — 1 and the area law for d > 1 for a large 
class of fermionic and bosonic c?-dimensional many-body 
Hamiltonians with short-range interaction at T = can 
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be unified through S%(L) oc [L d ~ l - l]/(d - 1) = ln 2 _ d L 
(i.e., nonextensive, InL for d = 1, and L d_1 for d > 1, 
area law) [45|], which would correspond to a large class 
(not yet completely identified) of fully entangled quan- 
tum systems. For all these systems, one could expect 
that a value of q exists such that S q (L) oc L d (i.e., exten- 
sive). In this paper our conjecture is verified for d = 1 
quantum spin chains and d = 2 bosonic systems and is 



therefore promising also in higher dimensions. 
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